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ABSTRACT 

We present scaling relations between stellar-mass (Af*) and the size of galaxies at 0.3 < z < 
3 for half- (i?5o) and 90 percent-light (i?9o) radii, using a deep if -band selected catalogue 
taken with the Subaru Telescope and MOIRCS in the GOODS-North region. The logarithmic 
slope R oc Af^-^^'^-^ is independent of redshift in a wide mass range of A/* ^ 10^ — 10^^ 
Mq, irrespective of galaxy populations (star-forming, quiescent). The offset change is < 50 
percent. Provided that optical light in the rest frame traces the stellar mass of galaxies, the 
universal relation demonstrates that the stellar mass was built up in galaxies over their cosmic 
histories in a similar manner on average irrelevant to galaxy mass. The small offset in each 
stellar mass bin from the universal relation shows weak size evolution at a given mass. There is 
a moderate increase of 30-50 percent for i?5o and i?9o for less massive galaxies (Af* < 10^° 
Mq) from z ^ 3 to z ^ 1, while the sizes remains unchanged or slightly decrease towards 
z ~ 0.3. For massive galaxies (Af^,, > 10^^ M©), the evolution is ^ 70 — 80 % increase in i?go 
from z ~ 3 to z ^ 0.3, though that in i?5o is weaker. The evolution of compactness factor, 
i?5o/i?9o, which becomes smaller at lower redshift, is suggestive of minor merging effect in 
the outer envelope of massive galaxies. 

Key words: galaxies: evolution - galaxies; fundamental parameters - galaxies: high-redshift 
- infrared: galaxies. 



; 1 INTRODUCTION 

. The scale size of galaxies is one of the fundamental parameters 

■ to elucidate the history of galaxy formation and evolution. The 

■ change of size and stellar-mass relations over cosmic time would 
pose strong constraints on models of galaxy evolution. The obser- 
vational relations between galaxy size and stellar mass have been 
studied in the local universe, based on the Sloan Digital Sky Survey 
(Shen et al. 2003; Bemardi et al. 2011). Using rest-frame optical 
bands, which presumably trace the distribution of stellar mass in 
galaxies, many studies have investigated galaxy sizes at higher red- 
shift as a function of stellar mass for massive galaxies (A/* > 10^" 
Mq). For example, the relations for galaxies at 0.2 < z < 1 were 
studied for late-type galaxies (Barden et al. 2005) and for early- 
type galaxies (Mcintosh et al. 2005; Trujillo, Ferreras, & de la Rosa 
2011). Damjanov et al. (2009) and Mancini et al. (2010) gave size- 
mass relations of massive galaxies (A/* > lO^''-^ M©) at z ~ 1.5. 
Williams et al. (2010) studied the relation with large samples of 
galaxies at z < 2. For higher redshifts to z ~ 3, size-mass rela- 
tions have been obtained for galaxies with M, > W^" Mq (e.g., 
Franx et al. 2008; Nagy et al. 201 1 ; Cassata et al. 201 1). 

Many studies have corroborated that massive galaxies at high 
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redshifts were much smaller than local galaxies with comparable 
mass (e.g., Daddi et al. 2005; Trujillo et al. 2006, 2007; Toft et al. 
2007; Zirm et al. 2007; Cimatti et al. 2008; Buitrago et al. 2008; van 
Dokkum et al. 2008, 2009, 2010; Akiyama et al. 2008; Damjanov 
et al. 2009; Carrasco et al. 2010; Cassata et al. 2010; Szomoru et 
al. 2010; van der Wei 201 1; Cassata et al. 201 1). At a fixed stellar 
mass, spheroidal galaxies were significantly more compact at high 
redshift and evolved with rapid increase of the effective radius by a 
factor ~4 or even larger from z ~ 2 (e.g., Buitrago et al. 2008; Car- 
rasco et al. 2010) and by a factor ~2 from z ^ 1 (e.g., van der Wei 
et al. 2008; Trujillo et al. 2011). The finding of compact massive 
galaxies with a high stellar velocity dispersion also supports their 
existence (van Dokkum 2009; van de Sande 201 1). It contrasts with 
the absence of such compact massive galaxies in the local universe, 
though several candidates have been found at 2: ~ 0.5 (Stockton et 
al. 2010) and in the local universe (Trujillo et al. 2009; Valentinuzzi 
et al. 2010). The findings demonstrate that massive galaxies have 
increased their size dramatically since 2 ~ 3 in a different manner 
from the evolution of less massive galaxies. 

However, other studies have reached contradictory conclu- 
sions. There is significant disagreement between the results of dif- 
ferent studies. Barden et al. (2005) found weak or no evolution in 
the relation between stellar mass and effective disc size for galax- 
ies with Mt > 10^" Mq since z ~ 1. For early-type galaxies at 
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z ^ 1, Mcintosh et al. (2005) showed that luminosity-size and stel- 
lar mass size relations evolve in a manner that is consistent with a 
completely passive evolution of the red early-type galaxy popula- 
tion. It is also shown that not all high-redshift early-type galaxies 
were compact and underwent dramatic size evolution (e.g., Toft et 
al. 2007; Zirm et al. 2007; Saracco et al. 2009, 2011; Mancini et 
al. 2010; Stott et al, 2011). From the study of surface brightness in 
rest-frame V and z bands at 2; < 3, Ichikawa et al. (2010) gave 
another evidence for no conspicuous evolution in galaxy sizes. 

As many previous studies show, the size evolution of galax- 
ies are still controversial. Any systematic errors in the observation 
or analyses could bias results towards such a significant evolution 
(e.g., Mancini et al. 2010; Hopkins et a. 2009b; Bouwens et al. 
2004). The origin of the discrepancy could be ascribed to redshift 
effects that more distant galaxies look more compact due to the 
difficulty of measuring envelopes at low surface density. The light 
from the outer portion of high-redshift galaxies is apt to be hid- 
den in noise for low S/N observations. As the consequence, effec- 
tive radii and total luminosity (or stellar mass) would be underesti- 
mated. On the other hand, very deep observations (e.g., Szomoru et 
al. 2010; Cassata et al. 2010; Law et al. 2012) or stacking methods 
to enhance the faint envelope of galaxies (e.g., Zirm et al. 2007; van 
Dokkum et al. 2008, 2010; van der Wei et al. 2008) have claimed 
that it is not the case. 

How significantly have the sizes of less-massive normal galax- 
ies evolved from the early universe to the current epoch? In this 
context, we look into the evolution of stellar-mass and size scaling 
relations on the basis of half- and 90 percent-light encircles, focus- 
ing on less massive galaxies (M* < 10^^ M©) at 0.3 < 2; < 3, 
using a deep /sT-selected galaxy catalogue. We infer that the outer 
radius is more influenced by merging effect than star formation or 
central activity. In §2, we describe the catalogue we used, which is 
among the deepest in the K band to date. The depth is crucial for 
studying galaxies of low-surface brightness or galaxies which are 
dimmed due to the cosmological expansion at high redshift to mea- 
sure the radius at faint outskirt of galaxies. The data analysis and 
the result for the size and stellar-mass relation are detailed in §3 and 
§4. The results are discussed in §5. Ichikawa et al. (2010) studied 
the evolution of surface brightness of galaxies at 2; < 3 in rest- 
frame V and z bands with the same data as the present study. We 
will discuss the consistency of the present result with their study. 

Throughout this paper, we assume Qm = 0.3, Qa = 0.7, and 
Ho = 70 km s~^ Mpc~^. We use the AB magnitude system (Oke 
& Gunn 1983; Fukugita et al. 1996). 



2 DATA 

We use the K-bsmd selected catalogue of the MOIRCS Deep Sur- 
vey (MODS) in the GOODS-North region (Kajisawa et al. 2009, 
hereafter K09; Kajisawa et al. 2011, Kll), which are based on our 
imaging observations in JHKs bands with MOIRCS (Suzuki et 
al. 2008) and archived data. Four MOIRCS pointings cover 70 per- 
cent of the GOODS-North region (103 arcmin^, hereafter referred 
as 'wide' field). One of the four pointings, which includes HDF-N 
(Williams et al. 1996), is the ultra-deep field of MODS (28 arcmin^, 
'deep' field). As the accuracy of background subtraction was highly 
demanded for the study of faint end of galaxies, the background 
was scrutinized and carefully subtracted (see Kll for the details). 
The surface brightness limit was extensively examined in Kl 1 . The 
t5?pical Icr siuface brightness fluctuations in one arcsec diameter 
were found to be ~27 mag arcsec"^ and ~27.5 mag arcsec"^ in 



K band for the wide and deep fields respectively. These are 1 ~ 2 
mag deeper than those of previous studies in K band. The depth is 
crucial for the study of low surface brightness features at the out- 
skirts of galaxies because of the strong dependence of cosmological 
dimming of surface brightness on redshift. 

For the total magnitude rriK vciK band, we use MAG_AUTO 
obtained by SExfi-actor (Bertin & Amouts 1996). The 85 ~ 90- 
percent completeness of the catalogue is rriK ~ 25 in the wide field 
and ~26 mag in the deep field. We exclude fainter galaxies from 
the present analysis. The FWHMs of the final stacked images are 
0.46 arcsec for the deep image and 0.53-0.60 arcsec for the wide 
image. The numbers of galaxies are 3555 and 6063, respectively. 

To obtain the stellar mass of MODS samples, K09 performed 
SED fitting of the multiband photometry (UBVizJHK, 3.6 fim, 
4.5 jjm, and 5.8 /im) with population synthesis models. We adopt 
the results with GALAXEV templates (Bruzual & Chariot 2003) 
and the Salpeter initial mass function (see K09 for more details). 
The stellar masses (M*) are obtained from the best-fit stellar mass- 
to-luminosity ratio in K band and scaled with the K-band flux. 
A detailed comparison between different mass estimators is given 
in K09. In the present analysis, the near-infrared data (3.6yum) of 
Spitzer/IRAC are available for most of the sample galaxies (96 per- 
cent), so that SED fitting is reasonably reliable for the photometry 
at rest-V^ (Acff — 0.55/im) out to redshift ~3. We assume that the 
size of the stellar system in galaxies is represented by the size mea- 
sured on -band images, which are the rest-frame optical to near- 
infrared wavelengths at z < 3. The difference of the sizes between 
optical and infrared light is discussed in §4. 

In the present catalogue, 209 galaxies are identified as X-ray 
sources. Among them, 61 sources are massive galaxies emitting 
hard X-rays with M* > 10^°^ M© at 2 < 2: < 4 (Yamada et 
al. 2009). One would expect galaxies with AGN emission to have 
smaller half-light radius. Considering possible effects on the size 
and mass estimates, we discard X-ray sources in what follows. 



3 ANALYSES 
3.1 Galaxy sizes 

As the definition of galaxy size, most previous studies have used the 
scale length (re), obtained by fitting the Sersic profile (Sersic 1968) 
to galaxy images. The accuracy of the scale length strongly depends 
on image quality and depth of observations. For bright galaxies, re- 
liable galaxy sizes can be measured from ground-based data (Tru- 
jillo et al. 2006; Franx et al. 2008; Williams et al. 2010). On the 
other hand, it is crucial to take an accurate measure of the profile for 
small objects near the resolution limit. In fact, Konishi et al. (201 1) 
applied two dimensional fitting for the present sample to obtain 
Sersic indexes with a single component for a morphological study 
of MODS galaxies. However, they found that reliable fitting was 
successful only for comparatively massive galaxies (M, > 10^'^ 
Mq) at 2 ~ 1 and very bright galaxies (K < 22.5) at higher 
redshift. The resolution of FWHM = 0.5 ~ 0.6 arcsec like the 
present sample would not allow reliable fitting with Sersic profile 
to fainter galaxies. Sersic parameters are sometimes degenerated 
between n and re, which depends on the surface-brightness limit 
of the images (e.g., Mancini et al. 2010; Stott et al. 2011), so that 
the reliable fitting would be difficult for faint galaxies. Moreover 
high-redshift galaxies exhibit a wide range of disturbed morpholo- 
gies (e.g., Kajisawa & Yamada 2001). Therefore, the fitting for less 
massive galaxies at high redshift could also be strongly influenced 
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by their more amorphous features in the shape. In addition, we must 
take into account that we are observing a mix of galaxy morpholo- 
gies. 

As such, we define the half-hght radius (rso) as the radius of 
a circular aperture, which encircles half the Tf-band light emitted 
from galaxies. As we here focus not just on massive (or bright) 
galaxies, but also on less massive galaxies up to 2; ~ 3, we prefer 
rso to Te for the study of the structural parameters. If a galaxy pro- 
file extends to infinity, encloses half the flux in the Sersic profile. 
However, as galaxies are supposed to have an edge at several times 
the scale length, does not always mean half-light radius. There- 
fore, in general, rso is smaller than re- In the same way, we de- 
fine 90 percent-light radius rgo. The 90-percent radius will give us 
more information at outer region of galaxies, where the size would 
be more influenced by merging effect, rso and rgo are obtained by 
SExtractor with PHOT_FLUXFRAC=0.5 and 0.9 

Figure [T] shows rso versus tuk for all objects in MODS. 
MODS catalogue lists many spectroscopically confirmed stars, 
which are also plotted in Fig. [T] As few spectroscopic data were 
available for stars fainter than rriK ~ 23.5 in the present region, 
we examined the reliability of rso using artificial stars. The arte- 
facts were convolved with a Moffat point spread function (/? = 3 
and FWHM=0.5 arcsec in the deep field and 0.6 arcsec in the 
wide field), and then buried in images as noisy as those of the 
deep and wide fields. SExtractor was applied to the artefacts in the 
same manner as for the MODS catalogue. The results are plotted in 
Fig. [T] The location of stars gives a clear boundary of unresolved 
galaxies. We define galaxies smaller than rso = 0.41 in the wide 
field rso = 0.30 in the deep field as unresolved galaxies. 

3.2 Monte Carlo study for the sizes with mock galaxies 

We first study the reliability of the galaxy size. Seeing and back- 
ground noise strongly affect the observed rso and rgo, especially 
for faint or small galaxies. To examine the effect on the size esti- 
mate, we generated mock galaxies with a 1/4-law or exponential 
profile extending to 3 times scale length. The images are then con- 
volved with a Moffat point spread function (/3 = 3 and FWHM=0.5 
arcsec in the deep field and 0.6 arcsec in the wide field), and buried 
in the simulated noise image. The galaxies were randomly gener- 
ated with various magnitudes and effective radii or scale lengths, 
and analyzed with SExtractor in the same manner as for the MODS 
galaxies. The results are shown in Fig.|2] Seeing and noise seriously 
change the observed radii specifically for galaxies with a 1/4-law 
profile. The observed rso is significantly smaller than the intrinsic 
values in general. The effect is stronger for fainter 1/4-law galaxies, 
whereas the effect is much smaller for exponential galaxies. The ap- 
parent smaller size than the intrinsic size would be due to the fact 
that the faint and large envelope of 1/4-law galaxies is hidden in 
background noise. We should be aware of this effect if the obser- 
vation is not deep enough for early-type galaxies (see also Hopkins 
et al. 2009b; Mancini et al. 2010). The effect will be discussed in 
more detail in the following section. It is noted that Williams et al. 
(2010) claimed no systematic effects on for rriK ^ 24, using the 
shallower MODS images. 

The radii, rso and rgo, for small galaxies are affected by 
smearing due to seeing effect and the effect depends on galaxy mor- 
phology and size. Therefore, we obtained Eq. 1 to correct the size 
effect by fitting the deviation of the size as shown in Fig. [2] (dash 
line). However, it is hard to examine the shapes of the present small 
or high-2 galaxies, so that we do not correct the effect which de- 
pends on morphology. 



rintrinsic — •\/^obscrvGd (At^. (1) 

We adopt Ar = 0.25 and 0.35 arcsec for rso and rgo, respectively. 

The size errors for mock galaxies are depicted in Fig. [5] as a 
function of observed magnitude. We should take into account that 
much smaller values are obtained for the size of 1/4-law galax- 
ies. Figure|4]is the difference between intrinsic and observed mag- 
nitudes. At the faint limit, observed magnitudes are fainter by 
1 ^ 1.5 mag for 1/4-law galaxies, while the effect is much smaller 
for exponential galaxies. These results suggest that shallower ob- 
servations tend to more underestimate size and flux (therefore, stel- 
lar mass). We discuss later the influence on the results due to the 
systematic errors. 

3.3 Size and stellar mass relations 

Using photometric (or spectroscopic, if available) redshift data, rso 
and rgo in arcsec are converted to the physical size in kpc (i?5o 
and -Rgo). The photo-z accuracy of the present sample is reason- 
ably good, Sz/{1 + z) = -0.011 (o- = 0.078) (KU), so that we 
expect no serious influence on our result due to the photo-z error 
We will confirm this later using galaxies with spectroscopic red- 
shifts. The results for _Rso and _Rgo are shown in Figs. [5] and |6] as 
a function of the stellar mass of galaxies for different redshift bins. 
The unresolved galaxies are displayed in another figure (Fig. O 
for reference. The results of effective radius (r^) for local galax- 
ies by Shen et al. (2003) are also shown in Fig.[5]for comparison. 
Size-mass relations for star-forming galaxies are consistent with 
that of late-type galaxies by Shen et al. (2003), although our defini- 
tion of half-light radius is different from that of Shen et al. (2003). 
It should be noted that it would be difficult to discuss the consis- 
tency for massive early-type galaxies because of the small number 
of such galaxies in our sample. 

Since our image quality is not high enough for classifying the 
galaxies into morphological classes at high redshifts, we divided 
the samples into quiescent and star-forming galaxy groups using a 
two-color diagnostic plot of rest-frame U -■ V and V ~ J colors 
(Williams et al. 2009). The adopted selection criteria for quiescent 
galaxies are as follows: 

iU-V)> 0.88{V - J) + c, (2) 

where U — V and V — J in the rest frame were obtained with the 
SED model fit to galaxies. The offset, c, is 0.69, 0.59, and 0.49 for 
0.0 < z < 0.5, 0.5 < 2: < 1.0, and 1.0 < 2; < 2.5, respectively. 
For z > 2.5 galaxies, we applied the offset for 1 < 2: < 2.5. 
Therefore the selection would be less reliable (see Williams et al. 
2009), though the number of such galaxies is very small. Additional 
criteria of (7 — V > 1.3 and 1^ — J < 1.6 are required for quiescent 
galaxies at all redshifts to prevent contamination from unobscured 
and dusty star-forming galaxies, respectively. 

3.4 Stellar-mass surface density 

The small differences in the slope and offset of the regression lines 
in different redshift bins in Figs. [5] and |6] would imply a univer- 
sal relation between the stellar mass and size of galaxies, irrele- 
vant to redshift. In this context, we plot the stellar-mass surface 
density (SMSD), which is defined as ^, = log(Af»/7ri?'^), in a 
single redshift bin as a function of stellar mass for all galaxies at 
0.25 < 2 < 3 (Figs.|9land|10t. The figures are essentially equiva- 
lent to the size-mass distribution. However, we can easily identify 
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Figure 1. Half-light radius (rso) and total magnitude (m/^:) for the MODS catalogue. The dots (red) and crosses (blue) shows the samples in the deep and wide 
fields, respectively. The vertical solid and dashed lines indicate the magnitude limit for each field. Filled (blue) squares and circles (red) are spectroscopically 
confirmed stars. Few stars fainter than 23.5 mag with spectroscopic data are available, so that we plot the results of SExtractor detection of artificial stars 
embedded in noise images (open blue squares for the wide and open red circles for the deep). The horizontal (blue) dash and (red) solid lines delineates 
unresolved sources. 



very compact or low-surface brightness galaxies in the figures, if 
they exist. For example, galaxies as compact (or diffuse) as 3 times 
smaller (larger) in size than the average should be located at ~ 1.0 
dex above (below) the average. Such low surface brightness galax- 
ies can be recognized at i\f* < 10^° M© in Fig. 1101 

In contrast, massive and compact galaxies are rare in the 
present sample. It is noted that most of massive galaxies with 
M* > 5 X 10^" M© should be resolved in all redshift bins of 
the deep field. We conjecture that very compact galaxies coalesce 
in the unresolved galaxies in Fig.l?] 



4 RESULTS 

4.1 Linear regression of size-stellar mass relation 

Since the strong correlations of 7?5o and Rqq with Af* in Figs. |5] 
and|6]are suggested, we obtain the least square fit between the size 
and mass of the galaxies with a linear regression, 

logii = a^(logM, - M) + logi?", (3) 



where r is 50 for half light radius or 90 for 90% light radius and 
is the radius at Kh = 10^^ Mq. M is 10 for all and star- 
forming galaxies, while AI — 11 for quiescent galaxies because 
the quiescent galaxies with < 10^" Mq were not observed at 
2 > 1.5 and their is statistically less robust due to the small 
numbers of the sample in narrow mass ranges at high redshits. We 
note that quiescent galaxies are located on average below the re- 
gression lines in all redshift bins, though the offset is smaller for 
Rgo than for 7?5o. The fact suggests that quiescent galaxies are 
more compact than star-forming galaxies at a give mass. The best- 
fit slopes (aso and ago) and offsets (R^q and i?go) of the regression 
analysis with mean errors for all, quiescent galaxies (QSG), and 
star-galaxies (SFG) are summarized in Table 1. The dispersion, a, 
of the linear fit is listed in the sixth and last columns. As the exclu- 
sion of unresolved sources could biases the results towards larger 
radii, we obtained in Table 2 the least-squares fit with the unre- 
solved galaxies in Fig|7] 

We depict the evolution of the slope and offset as a function 
of redshift in Fig.jS] While the slope for the quiescent galaxies are 
a little steeper than those for star-forming galaxies, the figure in- 
dicates that the slopes remain within ~ 0.1 — 0.2 and offsets do 
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Figure 2. (a) Half-light (rso) and (b) 90 percent-light (rgo) radii of mock galaxies embedded in the noise images. The abscissa is the intrinsic values and the 
ordinate is those observed with SExtractor Triangles (blue) and circles (red) indicate the model galaxies with exponential and 1 / 4-law profiles, respectively. 
Open symbols are the sources (rriK ^ 25) in the noise image of the wide field, while filled ones are those {rriK ^ 26) in the deep field. The sample galaxies 
larger than rso = 0.41 arcsec for the wide field and rso = 0-3 arcsec for the deep field are plotted above horizontal dot lines in (a). The dashed curves show 
the effect of point spread function Ar = 0.25 and 0.35 arcsec for rso and rgo, respectively. 
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Figure 3. Ratio of observed size to intrinsic for (a) rso and (b) rgo as a function of the observed magnitude for mock galaxies. The symbols are the same as 
in Fig. [2] 



not significantly change from z ~ 3 to 2: '-^ 0.3, irrespective of the 
sample selection except the quiescent population at z > 2 where 
the statistical error is very large. 

For SMSD in Figs.|9]and[T0l we define the regression as 

^1, = a^log(M. - 10) + Atr°, (4) 

where is SMSD at M. = 10^° Mq. The results are shown in 
Table 3. 



4.2 Error estimate 

The errors in stellar mass mainly originated from SED fitting and 
photometric errors of the observations. In the course of fit- 
ting to obtain the best SED model with various parameters (e.g., 
star-formation time scale, photometric redshift for galaxies with 
no spectroscopic redshift available, age, extinction, and metallic- 
ity), the probability distributions of stellar mass can be calculated, 
where the photometry and photometric-redshift errors are included. 
The error for the observed SMSD is dominated by the error for the 
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Figure 4. Difference between the intrinsic and observed magnitudes (observed minus intrinsic) as a function of the observed magnitude. The symbols are the 
same as in Fig.|2] 



stellar mass. See more details for the error estimate in K09 and 
Ichikawa et al. (2010) 

The size and magnitude errors in accordance with galaxy mag- 
nitude have been examined using mock galaxies (Figs.[3]and|4l(. If 
all galaxies have a shape of 1/4-law, the error would be significant. 
Recalling the size error due to the morphology strongly depends on 
magnitude, we examine the effect by confining the sample to bright 
galaxies, where the systematic error becomes smaller. We selected 
the samples with the magnitude limits, mx = 24 (23), 25(24), 
26(25) for the deep (wide) fields, then obtained the regression lines 
again. The results are compared in Table 3. It should be noted that 
the section does not change the result. 

One would be concerned about the reliance on photometric 
redshift estimates for high redshift galaxies (e.g., Mosleh et al. 
2011). In order to examine the error, we selected only the galax- 
ies with spectroscopic redshift available and compared the result 
with the photo-2: sample. Although the spectroscopic samples are 
limited to lower redshift, the result does not change our conclusion 
(Table 3). 

Due to possible color gradients in galaxies, it would be best 
to use images of the same rest-frame band with comparable S/N 
for measuring the size of galaxies at all redshifts. Nevertheless, we 
used the radius of galaxies measured on K-hanA image, because 
it is deepest among the images we used and because it is the rest- 
frame optical or longer band at z < 3. If the central region of galax- 
ies are younger (bluer) than the outer region, the mass-weighted 
radius could be larger than the luminosity-weighted radius. If it is 
older (redder), the result would be vice versa. Star forming galaxies 
sometimes show strong morphological variation between observed 
wavelengths. Bond, Gawiser & Koekemoer (2011) reported that 



this was not generically accompanied by a large difference in half- 
light radius. Barden et al. (2005) measured the disk scale lengths of 
local galaxies in various bands. The average size is about 10 percent 
larger in V band than in K band. On the other hand, MacArthur 
& Courteau (2003) showed the contrary result that the distribu- 
tion of disc scale lengths was a decreasing trend with increasing 
wavelengths (see also Cassata et al. 2010 for early-type galaxies at 
z ~ 2). 

To investigate the effect of the color gradient, we compare the 
90 present-light radius of galaxies measured on the nearest rest- 
frame V band for each redshift bin. The sizes of the present sam- 
ples were obtained on ACS J, z, and MOIRCS J, H bands and 
compared with that of K band in Fig. [TT] We used ACS images 
binned to 0.117 arcsec per pixel, keeping the original image reso- 
lution. Fig. 1 1 1 [demonstrates that the sizes of ACS in I and z bands 
are systematically ~ 25 percent smaller than that in K band. 

The convolution of the ACS J-band image with PSF and see- 
ing enhances low surface brightness details. If the ACS images are 
convolved with a gaussian (FWHM=2 pixel, 0.234 arcsec), the dif- 
ference is decreased to ~15 percent. As the convolution enhances 
the galaxy edge of low surface brightness, it tends to give a larger 
galaxy size (see also fig. 7 of Mancini et al. 2010). It should be 
noted that the J and H images gives larger image size for smaller 
galaxies. The size strongly depends on the depth and the seeing 
size. 

4.3 Size evolution 

Although we have found the universal relation between Af, and 
SMSD (or size), which does not depend strongly on redshift, we 
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log{M.) (Mq) 

Figure 5. Physical half-light radius (-R50), corrected for PSF effects (Eq.l), as a function of the stellar mass (M*) of the sample galaxies. The galaxies are 
classified into quiescent (red) and star-forming galaxies (blue) (see text) and the fields (wide and deep) as noted in the bottom left frame. Solid line is the 
linear regression obtained for all samples in each redshift bin. Red dash-dot and blue dot lines are the regressions for quiescent and star-forming galaxies, 
respectively. Black dashed and dash-dot lines in the 0.25 < z < 0.5 bin are the results of effective radius (re) for local late- and early-type galaxies (Shen et 
al. 2003), respectively. The horizontal dot lines are the size limits of the present observation for the wide (upper) and deep (lower) fields at the lower redshift 
boundary for each redshift bin. 



have good reason to expect offsets from the relation for some 
galaxy populations (e.g., massive compact galaxies) in certain 
mass and redshift ranges. The regression analysis could be more 
weighted on the numerous less massive populations. Therefore, it 
is likely that there are small (but not highly significant) offsets be- 
tween galaxy populations. The small offsets would account for the 
size evolution of such populations. Using the M^-R relation ob- 
tained for the galaxies at 0.25 < z < 0.5 (Tables 1 and 2) as a 
reference, we examined the deviation of the median size of galax- 
ies in the mass and redshift bins from the reference. The evolution 
of the massive (M* > lO^^^'^M©) quiescent galaxies, defined in a 
U — V and V — J diagram (Williams et al. 2010), is also obtained. 
(We use median values to avoid unreasonable contributions from 
outliers with large deviation.) We show the results in Fig.[T2l where 



unresolved galaxies are included. In addition, the median values of 
R50/R90, which represent a sort of compactness of galaxies, are 
depicted in the figure to see the evolution of the compactness as a 
function of redshift. We note that the result with resolved galaxies 
are in good agreement with that with unresolved galaxies. 



5 DISCUSSION 

Using the deepest i('-band image, we obtained scaling relations 
between stellar-mass (A/,) and size (or equivalently SMSD) as 
a function of redshift for galaxies at 0.3 < z < 3 in a wide 
mass range of Ah ^ 10* — lO'^^ M©. We defined the radii en- 
circling half- (i?5o) and 90 percent-light (Rgo) in a circular aper- 




log(M.) (Mq) 

Figure 6. Same as Fig. [5] but for 90 percent-light radius (Rgo)- 



ture. The depth was found to be crucial for proper measure of the 
size and mass for less massive galaxies dimmed due to cosmic ex- 
pansion. We conclude that there is no strong evidence for the size 
evolution at a given mass over the redshift range, irrespective of 
galaxy populations (star-forming, quiescent). The size-mass rela- 
tion (R oc A/") have a universal slope (a ~ 0.1 — 0.2) and small 
offsets (< 50 percent). In other words, as galaxies grow in mass 
through star formation or merging process, their sizes evolve in 
such a way that galaxies in general move along the scaling relation. 
The universal relation demonstrates that the stellar mass in galax- 
ies with the same stellar mass was built up on average in a similar 
manner over cosmic times. The trend is insensitive to the stellar 
populations. Our result is comparable with the oc M°'^^ found 
by Shen et al. (2003) for low-mass galaxies {M, < 10^°-^ Mq). 
On the other hand, it is in contrast to the steeper relation for massive 
galaxies obtained by previous studies (e.g., Shen 2003; Bemardi et 
al. 2011). The scaling relations are much tighter (a ~ 0.12-0.17 
dex) than those (cr ^ 0.3-0.5 dex) of previous studies in the local 



universe (e.g., Shen et al. 2003) or at high redshifts (e.g., Franx et 
al. 2008). 

The weak growth of the galaxy size, irrelevant to galaxy mass, 
is in disagreement with the scenario that more massive galaxies 
rapidly changed their sizes. If massive high- 2; galaxies are sev- 
eral times smaller than local galaxies with comparable mass, they 
should be located well above the stellar mass vs. SMSD relations. 
Several candidates of such compact galaxies with A/, > 10^^ pos- 
sibly coalesce as unresolved galaxies in Figs. |9] and [TO] However, 
such galaxies are found to be few. Our finding is in contrast with 
the previous results for quiescent early-type massive galaxies (e.g., 
Zirm et al. 2007; Williams et al. 2010). With regard to SMSD, the 
average SMSDs in R50 for our sample galaxies with A/, > 10^" 
M0 are found not to evolve with redshift (log fiR^g — 8.53 - 
8.70). The average SMSD for all galaxies at 0.25 < ^ < 3.0 is 
8.60 ± 0.01, which is very consistent with the result, 0.85 ± 0.03, 
for disk galaxies with Af, > obtained by Barden et al (2005), 
though our sample could be a mixture of disk and elliptical galax- 
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log(M.) (Mq) 

Figure 7. Same as Fig.|5] but for unresolved galaxies. The same regression lines for all sample in each redshift bin as in Fig.|5]are shown for reference. The 
dotted lines are the size limits for the wide (upper) and deep (lower) fields at the upper redshift boundary for each redshift bin. 



ies. However, due to small volume of the present study, the results 
could be subject to statistical uncertainties. 

The weak size evolution presented in the present study is 
reconciled with the surface brightness (SB) evolution. Ichikawa 
et al. (2010) presented a universal linear correlation between the 
stellar mass and SB in rest-frame V and z bands for galaxies at 
0.3 < z < 3, using the same sample as the present study. The 
correlation has a nearly constant slope, independent of redshift 
and color of galaxies in the rest-z frame. In contrast, SB shows 
a strong dependence on redshift for a given stellar mass. It evolves 
as (l + 2)-2-0~-0-«.The increase in the luminosity with redshift is 
estimated by using the expected luminosity evolution from a single 
burst at 2 = 4 for massive galaxies (A/, > 10^") and constant star 
formation for blue samples (U ~ V < 0) for less massive galax- 
ies. The redshift dependence of SB evolution is well explained by 
the pure luminosity evolution of galaxies out to z ~ 3 without the 
need of the size evolution, which supports the results for young 



early-type galaxies (e.g., Saracco et al. 201 1) and for disk galaxies 
at 2 < 1 (Barden et al. 2005) . 

Nevertheless, there are some populations having a small off- 
set from the universal scaling relation. We see more or less evolu- 
tion in the size and compactness for such populations as a func- 
tion of redshift in Fig. [T2l The figure suggests a moderate in- 
crease of 30-50 percent for R50 and Rgo for less massive galax- 
ies (Mt < 10^'' Mq) from 2 ~ 3 to z ~ 1, while the sizes re- 
mains unchanged or slightly decrease towards z ~ 0.3. For mas- 
sive galaxies (A/, > 10^^ M0), the evolution is ~ 70 — 80 percent 
in Rgo from 2 ~ 3 to 2 ~ 0.3, though that in Ri,o is weaker. Rqq 
evolved as oc (1 + 2)"° ■^ It is noted that Trujillo et al. (2006), 
using the ground-based K-bznd data comparable with our depth, 
but in a smaller region (6.3 arcmin^), concluded that there was no 
evidence for significant evolution for the size-mass relation and 
that there was small increase (29 percent) since 2 ~ 2.5 for the 
most massive bin M* > 4 x 10^° M©. In our sample, the average 
R50 for the most massive bin of M* = 10^^-10^^ Mq is larger 
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Table 1. Linear fit of radius to stellar mass for resolved galaxies. 



R50 R90 



Redshift 


Sample 


N 






(kpc) 


a 


N 




190 


(kpc) 


a 


0.25 ^ 


z 


< 3.00 


all 


5259 


0, 


,101 ± 0.002 


0.508 ±0.002 


0.130 


5236 


0, 


,103 ±0.003 


0.874 ± 0.002 


0.148 








QSG 


408 


0, 


,126 ± 0.009 


0.543 ±0.007 


0.110 


408 


0, 


,143 ±0.008 


0.955 ±0.006 


0.099 








SFG 


4851 


0, 


,113 ± 0.003 


0.520 ±0.002 


0.129 


4828 


0, 


,108 ±0.003 


0.880 ± 0.003 


0.151 


0.25 !S 


z 


< 0.50 


all 


543 


0, 


111 ± 0.006 


0.477 ± 0.009 


0.133 


542 


0, 


,114 ±0.008 


0.844 ±0.011 


0.168 








QSG 


28 


0, 


,129 ± 0.020 


0.537 ±0.026 


0.098 


28 


0, 


,154 ±0.021 


0.971 ±0.027 


0.100 








SFG 


515 


0, 


,115 ± 0.007 


0.486 ±0.010 


0.133 


514 


0, 


,114 ±0.009 


0.845 ±0.012 


0.170 


0.50 


z 


< 0.75 


all 


717 


0, 


,101 ± 0.006 


0.511 ±0.006 


0.119 


714 


0, 


,115 ±0.006 


0.883 ± 0.006 


0.124 








QSG 


70 


0, 


,133 ±0.018 


0.548 ± 0.019 


0.126 


70 


0, 


,149 ±0.015 


0.970 ±0.016 


0.107 








SFG 


647 


0, 


,119 ±0.007 


0.535 ± 0.007 


0.113 


644 


0, 


,122 ±0.007 


0.894 ± 0.008 


0.124 


0.75 ^ 


z 


< 1.00 


all 


1083 


0, 


,086 ± 0.005 


0.527 ± 0.004 


0.123 


1079 


0, 


,091 ± 0.006 


0.900 ± 0.005 


0.151 








QSG 


134 


0, 


,118 ± 0.017 


0.558 ±0.012 


0.104 


134 


0, 


,151 ±0.015 


0.983 ±0.011 


0.095 








SFG 


949 


0, 


,105 ± 0.006 


0.547 ± 0.005 


0.121 


945 


0, 


,098 ±0.007 


0.908 ± 0.006 


0.156 


1.00 


z 


< 1.25 


all 


888 


0, 


,097 ± 0.006 


0.531 ±0.005 


0.130 


885 


0, 


,100 ±0.006 


0.895 ±0.005 


0.124 








QSG 


83 


0, 


,153 ± 0.024 


0.553 ±0.014 


0.105 


83 


0, 


,176 ±0.019 


0.960 ±0.011 


0.085 








SFG 


805 


0, 


,122 ± 0.007 


0.554 ±0.006 


0.128 


802 


0, 


,114 ±0.007 


0.909 ± 0.006 


0.125 


1.25 


z 


< 1.50 


all 


519 


0, 


,100 ±0.008 


0.538 ±0.006 


0.127 


514 


0, 


,082 ±0.011 


0.903 ± 0.008 


0.161 








QSG 


32 


0, 


,125 ±0.037 


0.570 ± 0.020 


0.089 


32 


0, 


,123 ±0.033 


0.958 ±0.018 


0.079 








SFG 


487 


0, 


,114 ±0.009 


0.550 ± 0.007 


0.127 


482 


0, 


,089 ±0.012 


0.909 ± 0.009 


0.165 


1.50 ^ 


z 


< 2.00 


all 


555 


0, 


,123 ±0.010 


0.514 ±0.006 


0.128 


554 


0, 


,120 ±0.010 


0.869 ± 0.006 


0.128 








QSG 


28 


0, 


,166 ± 0.074 


0.518 ±0.027 


0.121 


28 


0, 


,158 ±0.058 


0.916 ±0.021 


0.094 








SFG 


527 


0, 


,148 ± 0.010 


0.528 ±0.006 


0.124 


526 


0, 


,136 ±0.011 


0.878 ± 0.006 


0.128 


2.00 < 


z 


< 2.50 


all 


700 


0, 


,128 ± 0.009 


0.468 ± 0.005 


0.121 


696 


0, 


,124 ±0.010 


0.829 ±0.006 


0.146 








QSG 


27 


0, 


,250 ±0.083 


0.484 ±0.021 


0.109 


27 


0, 


,233 ± 0.074 


0.881 ±0.019 


0.097 








SFG 


673 


0, 


,146 ± 0.009 


0.473 ± 0.005 


0.119 


669 


0, 


,135 ±0.011 


0.832 ± 0.006 


0.147 


2.50 ^ 


z 


< 3.00 


all 


254 


0, 


,127 ±0.017 


0.438 ± 0.008 


0.130 


252 


0, 


,099 ±0.022 


0.818 ±0.011 


0.166 








QSG 


6 


0, 


,480 ± 0.070 


0.479 ±0.019 


0.040 


6 


0, 


,416 ±0.062 


0.855 ±0.017 


0.035 








SFG 


248 


0, 


,140 ±0.018 


0.442 ± 0.009 


0.129 


246 


0, 


,109 ±0.024 


0.822 ±0.011 


0.167 



by ~ 20 — 30 percent at z ~ 0.3 than at z ^ 2.5, which is in 
good agreement with the result of Trujillo et al. (2006), whereas it 
is much weaker than those of other studies for galaxies of similar 
mass (e.g., Franx et al. 2008; Buitrago et al. 2008; van der Wei et 
al. 2009; Williams et al. 2010). Using new deep WFC3 data, Cas- 
sata et al. (2011) and Law et al. (2012) have given new evidences 
of compact galaxies at high-z. Although the data in H band is shal- 
lower by 1 mag arcsec^^ than our A'-band data, their conclusion 
would be robust. It would be possible that the compact galaxies are 
mingled in unresolved galaxies of the present sample or that the 
difference of the definition for half-right radius {r^ and rso) gives 
the contradict result. 

The ratio, R50/R90, in Fig. 1121 would give a clue to under- 
standing the evolution of the compactness of galaxies. While mid- 
dle and low mass galaxies shows no evolution in the compactness, 
the galaxies in the most massive bin are more compact at higher 
redshift. Recalling weaker evolution in R50 than in Rgo, we infer 
that the compactness evolution is ascribed mainly to the expansion 
of the outer rim of massive galaxies. In other words, at a given 
mass, massive galaxies in the local universe are more influenced 
by mergers or star formation at the outer envelopes than those at 
high redshifts. It would be worth noting that the galaxies in the 
massive bin of the deep field are complete in sampling over the 
present redshift range and their R^o and iigo are well resolved in 
the present PSF. However, the present observation samples a com- 
paratively small volume and therefore, the results could be subject 
to statistical uncertainties for massive galaxies of low number den- 
sity due to field variances. 



The size evolution is often used to advocate the merging pro- 
cesses in the hierarchical paradigm of galaxy formation and evo- 
lution. Minor dry merger is a plausible mechanism for weak size 
evolution (e.g., Guo & White 2008; Bezanson et al. 2009; Naab, 
Johansson, & Ostriker 2009; Hopkins et al. 2010). In that the size 
evolution is stronger in massive galaxies, our findings are quali- 
tatively in agreement with previous studies based on simulations 
(e.g., Boylan-Kolchin, Ma, & Quataert 2006; Hopkins et al. 2009a). 
Naab et al. (2009) showed that minor mergers or the accretion of 
relatively low-mass satellites may be the main driver for the late 
evolution of sizes of massive early-type galaxies. Somerville et 
al. (2008) showed based on a CDM model of disc formation with 

> 10^" Mq that the average size of discs at fixed stellar mass 
was about 50 percent larger than that at 2 ~ 3. The predicted evo- 
lution in the mean size at a fixed stellar mass since 2: ~ 1 is about 
15 ~ 20 percent, which is comparable with our observations. 

The newly accreted small galaxies preferentially populate the 
outer region of massive galaxies. The small size growth of the 
present result is plausibly accounted for by minor mergers or the 
accretion of relatively low-mass satellites (e.g., Naab et al. 2009; 
Hopkins et al. 2009a). Naab et al. (2009) and Fan et al. (2010) 
showed the fractional variation of the gravitational radius of the 
main galaxy after A'^ minor-merger events with R oc as 

where r f and are final and initial galaxy sizes, ry is the fractional 
ratio of merging galaxy to the main galaxy. It takes ~1 (~3) minor- 



A universal stellar-mass and size relation 1 1 



Table 2. Linear fit of radius to stellar mass with unresolved galaxies 



R50 -R90 

Redshift Sample N aso log^so <^ ^ «90 log-^go 

(kpc) (kpc) 



0.25 ^ 2 < 3.00 


all 


6532 


0, 


.131 ±0.003 


0.484 ± 0.002 


0.161 


6502 


0, 


,125 ±0.003 


0.852 ± 0.002 


0.161 




QSG 


445 


0, 


.132 ±0.008 


0.536 ± 0.007 


0.116 


445 


0, 


,147 ±0.008 


0.946 ± 0.006 


0.108 




SFG 


6087 


0, 


.144 ±0.003 


0.497 ± 0.003 


0.162 


6057 


0, 


,130 ±0.003 


0.859 ± 0.003 


0.164 


0.25 ^ 2 < 0.50 


all 


646 


0, 


,136 ± 0.007 


0.475 ±0.010 


0.156 


645 


0, 


,133 ±0.008 


0.838 ±0.012 


0.181 




QSG 


31 


0, 


,165 ± 0.023 


0.554 ±0.034 


0.129 


31 


0, 


,180 ±0.022 


0.983 ± 0.033 


0.124 




SFG 


615 


0, 


,141 ± 0.007 


0.485 ±0.011 


0.156 


614 


0, 


,133 ±0.008 


0.841 ±0.013 


0.183 


0.50 ^ 2 < 0.75 


all 


858 


0, 


,128 ± 0.006 


0.503 ±0.007 


0.136 


855 


0, 


,137 ±0.006 


0.875 ±0.007 


0.134 




QSG 


73 


0, 


,137 ±0.017 


0.546 ± 0.019 


0.127 


73 


0, 


,155 ±0.015 


0.968 ± 0.016 


0.111 




SFG 


785 


0. 


,150 ±0.007 


0.532 ± 0.008 


0.133 


782 


0, 


,149 ±0.007 


0.890 ± 0.008 


0.135 


0.75 ^ 2 < 1.00 


all 


1313 


0, 


,129 ±0.006 


0.512 ±0.005 


0.155 


1308 


0, 


,125 ±0.006 


0.886 ± 0.005 


0.161 




QSG 


141 


0, 


,119 ±0.015 


0.557 ±0.012 


0.102 


141 


0, 


,158 ±0.014 


0.984 ±0.011 


0.094 




SFG 


1172 


0, 


,152 ± 0.007 


0.5.36 ±0.006 


0.157 


1167 


0, 


,135 ±0.007 


0.898 ±0.007 


0.166 


1.00 2 < 1.25 


all 


1082 


0, 


,142 ± 0.007 


0.516 ±0.006 


0.164 


1077 


0, 


,129 ±0.006 


0.881 ±0.006 


0.146 




QSG 


88 


0, 


,150 ±0.022 


0.550 ±0.013 


0.103 


88 


0, 


,171 ±0.018 


0.958 ±0.011 


0.084 




SFG 


994 


0, 


,174 ±0.008 


0.546 ± 0.007 


0.164 


989 


0, 


,146 ±0.007 


0.898 ±0.007 


0.148 


1.25 ^ 2 < 1.50 


all 


675 


0, 


,145 ± 0.010 


0.515 ±0.008 


0.159 


668 


0, 


,120 ±0.011 


0.879 ± 0.009 


0.176 




QSG 


35 


0, 


,152 ±0.045 


0.562 ± 0.025 


0.111 


35 


0, 


,141 ±0.037 


0.951 ± 0.020 


0.091 




SFG 


640 


0. 


,163 ±0.011 


0.530 ± 0.009 


0.160 


633 


0. 


,130 ±0.012 


0.887 ± 0.010 


0.179 


1.50 ^ 2 < 2.00 


all 


685 


0, 


,166 ±0.012 


0.487 ± 0.007 


0.174 


684 


0, 


,133 ±0.009 


0.847 ± 0.006 


0.139 




QSG 


36 


0, 


,077 ± 0.069 


0.475 ± 0.024 


0.128 


36 


0, 


,044 ± 0.064 


0.865 ± 0.023 


0.120 




SFG 


649 


0, 


,206 ±0.013 


0.509 ± 0.008 


0.171 


648 


0, 


,159 ±0.010 


0.860 ± 0.006 


0.138 


2.00 ^ 2 < 2.50 


all 


923 


0, 


,187 ±0.009 


0.428 ± 0.005 


0.152 


918 


0, 


,165 ±0.010 


0.796 ± 0.005 


0.156 




QSG 


32 


0, 


,281 ± 0.068 


0.477 ±0.019 


0.103 


32 


0, 


,284 ± 0.062 


0.870 ±0.017 


0.094 




SFG 


891 


0, 


,208 ± 0.010 


0.436 ±0.005 


0.151 


886 


0, 


,179 ±0.011 


0.801 ±0.005 


0.157 


2.50 ^ 2 < 3.00 


all 


350 


0, 


,183 ± 0.017 


0.400 ± 0.009 


0.150 


347 


0, 


,139 ±0.020 


0.779 ±0.010 


0.175 




QSG 


9 


0, 


,251 ±0.133 


0.433 ± 0.039 


0.097 


9 


0, 


,157 ±0.148 


0.799 ± 0.043 


0.109 




SFG 


341 


0, 


,206 ±0.018 


0.409 ± 0.009 


0.149 


338 


0, 


,156 ±0.022 


0.786 ± 0.010 


0.176 



Table 3. Linear fit of surface mass density of galaxies in 0.25 < 2 < 3.0. 



Selection" 


Sample 




iir,u 














N 


M 
<0 


(Mq kpC-2) 


a 


N 


"90 


(Mq kpc-2) 


a 


TTiK =S 26(25) 


all 


5259 


0.80 ±0.00 


8.19 ± 0.00 


0.26 


5236 


0.79 ± 0.01 


7.71 ± 0.00 


0.30 




QSG 


408 


0.75 ±0.02 


8.37 ±0.01 


0.22 


408 


0.71 ± 0.02 


7.83 ± 0.01 


0.20 




SFG 


4851 


0.77 ±0.01 


8.16 ±0.00 


0.26 


4828 


0.78 ±0.01 


7.70 ±0.01 


0.30 


with unresolved galaxies 


all 


6532 


0.74 ± 0.01 


8.23 ± 0.00 


0.32 


6502 


0.75 ± 0.01 


7.75 ± 0.00 


0.32 




QSG 


445 


0.74 ± 0.02 


8.39 ± 0.01 


0.23 


445 


0.71 ± 0.02 


7.86 ± 0.01 


0.22 




SFG 


6087 


0.71 ±0.01 


8.21 ±0.01 


0.32 


6057 


0.74 ±0.01 


7.74 ± 0.01 


0.33 


rriK s£ 25(24) 




3980 


0.80 ±0.01 


8.19 ± 0.00 


0.25 


3972 


0.80 ± 0.01 


7.71 ± 0.00 


0.28 


rriK < 24(23) 




2378 


0.82 ±0.01 


8.20 ±0.01 


0.25 


2378 


0.79 ± 0.01 


7.71 ± 0.00 


0.23 


spectroscopic redshift 




1610 


0.81 ±0.01 


8.16 ±0.01 


0.24 


1610 


0.77 ±0.01 


7.69 ± 0.01 


0.22 



Selection for the deep field (parentheses are the sample for the wide field) 



merger events with a = 0.15 to increase the radius by ~ 20 (~ 70) 
percent for massive galaxies since z ^ 3, provided that the merger 
mass ratio is 1 : 10. It would not be unreasonable that massive galax- 
ies experienced such a small number of minor merging since z ~3 
(e.g., Btmdy et al. 2009; Lopez-Sanjuan et al. 2011). 

Deep observations will be important in constraining the ex- 
act amount (or lack thereof) and distribution of merging galaxies, 
and how galaxies built up with redshift. In this context, our find- 
ing demonstrates that minor mergers in massive system built up an 
envelope of lower siuface density materials. Deeper imaging ob- 



servations in a wider filed with high spatial resolution and consis- 
tent analyses from the local universe to high redshift will give con- 
straints on the compactness and the amount of low-surface bright- 
ness material at the outer envelope as a function of redshift to im- 
prove our understanding of galaxy formation. 
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Figure 8. Evolution of the slope and offset (Tables 2 and 3) as a function of redshift. Circles with solid line and crosses with dashed lines are the results for 
resolved galaxies and those with unresolved galaxies. Black, red, and blue symbols are the samples for all, quiescent, and star-forming galaxies, respectively. 
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Figure 10. Same as Fig.|9] but for 90 percent-light radius (-Rgo). 



Sersic J.-L. 1968, Atlas de Galaxias Australes (Cordoba: Obs. As- 
tron.) 

Somerville R. S., et al. 2008, ApJ, 672, 776 
Szomoru D., 2010, ApJ.,714, L244 

Stott J. P., Collins C. A., Burke C, Hamiltoii-Morris V., Smith 

G. P., 201 1, MNRAS, 414, 445 
Shen S., Mo H. J., White S. D. M., Blanton M. R., Kauffmanii G., 

Voges W., Brinkmaiin J., Csabai I., 2003, MNRAS, 343, 978 
Stockton A., Shih H-Y., Larson K., 2010, ApJ, 709, L58 
Suzuki R., et al. 2008, PAS J, 60, 1347 
Toft S., et al., 2007, ApJ, 671, 285 
Trujillo I., et al. 2006, ApJ, 650, 18 

Trujillo I., Conselice C. J., Bundy K., Cooper M. C, Eisenhardt 
P, EUis R. S. 2007, MNRAS, 382, 109 



Trujillo I., Cenarro A. J., de Lorenzo-Caceres A., Vazdekis A.„ de 

la Rosa I. G., Cava A., ApJ., 692, LI 18 
Trujillo I., Ferreras I., de la Rosa I. G., 201 1, MNRAS, 415, 3903 
Valentinuzzi T. et al. 2010, ApJ, 712, 1232 
van der Wei A., Bell E. R, van den Bosch F. C, Gallazzi A, Rix 

H-W., 2009, ApJ., 698, 1232 
van der Wei A., Holden B. P., Zirm A. W., Franx M., Rettura A., 

lUingworth G. D., Ford H. C, 2008, ApJ, 688, 48 
van der Wei A., et al., 201 1, ApJ, 730, 38 
van de Sande J., et al. 201 1, ApJ, 736, L9 
van Dokkum P G., et al., 2008, ApJ, 677, L5 
van Dokkum P G., Kriek M., Franx M., 2009, Nat, 460, 717 
van Dokkum P G., et al., 2010, ApJ, 709, 1018 
Williams R. E., et al. 1996, AJ, 112, 1335 



A universal stellar-mass and size relation 15 



o 

05 



o 



1 - 




2 

1 h 



1 F 


2 



1 ' 1 ' 1 ' 1 ' * 

/ 

0.25<z<0.50" 


1 1 ' 1 ' 1 ' 

z 

\- ■ 0.5<z<0.75 " 


— \ — i" ' V " i" "i — ^—^ — ^ — 


— 4 — i"- " \"' """i "'V" 1 — — 


z 

0.75<z<1.0 " 


J 

. .• 1.0<z<1.25 " 


> 1 ■ ■ 1 h — ^ ^ ^ ^ 


h — ^-H 1""' "l " i" '"1 ^ 


. ■ -. J 

' ' ■ 1.25<z<1.5 ' 

L- ■.::-;Kmi^^:^;.i:A . .. , J 


1.5<z<2.0 

;. . j<':>!;:Hi^ • •■■ , • 




— -■ — ''^'■y ^''b/i'iTy ■• — ■ 

■ • • •.v i v^m^^^ ■ . • • • 

'■•••■■'•.**•• '-i'"'!^: • ; 

-1,1,1,1, 


\ ^ ^ ^ ^ ^ 

~ J. 2.0<z<2.5 


0.5 1 1.5 2 


- • ■ - — ^ 

...... 

1 , ■ 1 , 1 , 1 



0.5 



1.5 



TqJK) (arcsec) 



90 
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Figure 12. Evolution of galaxy sizes as a function of redshift in stellar mass bins. The ordinate is the difference of the median sizes, (a) i?5o and (b) -R90, 
from the regression line obtained for the galaxies in 0.25 < z < 0.5 (Table 2). (c) Median size ratio R^o/Rgo. Um'esolved galaxies are included. The eiTor 
bar is the mean error of the average. 



